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SECTION  I 
INTRODUCTION 

Certain  re-entry  erosion  experiments  require  the  acceleration  of 
80-micron  size  particles  to  velocities  of  several  thousand  meters  per 
second.  A  suitable  way  to  produce  such  high  speed  particles  is  to 
accelerate  them  in  an  expanding  gas  flow.  The  fluid  dynamic  design  of 
the  necessary  expansion  nozzles  encounters  some  difficulties.  Eighty- 
micron  particles  are  large  as  far  as  particle  acceleration  in  a  gas 
flow  is  concerned  requiring  a  long  residence  time  and  thus,  a  long  path 
length  in  the  gas  for  the  acceleration  process.  Practical  path  lengths 
are  in  the  order  of  1  foot.  However,  to  obtain  a  sufficiently  high  gas 
velocity,  helium  at  high  pressure  and  temperature  is  used  as  an  operating 
medium.  The  scarcity  of  helium  and  the  great  extent  of  the  electric 
equipment  associated  with  the  heating  of  the  gas  (1  megawatt  arc  heater) 
hold  the  nozzle  mass  flow,  i.e.,  the  nozzle  diameter  down  (4  mm  throat 
diameter).  The  resulting  expansion  nozzles  are  very  slender  with  the 
consequence  that  the  expanding  gas  is  exposed  to  considerable  friction 
on  the  walls.  Determining  the  wall  friction  influence  becomes  a  major 
portion  of  the  design  process. 

For  the  present  nozzle  design  a  calculation  method  is  used  which 
expresses  the  expansion  process  in  a  given  duct  in  terms  of  the  occurring 
fluid  forces.  This  method  allows  a  ready  inclusion  of  wall  friction  forces. 
To  obtain  a  particularly  realistic  account  of  the  wall  friction  influence, 
this  method  is  applied  in  steps  to  consecutive  segments  of  the  expansion 
nozzle  with  each  segment  taken  small  enough  that  (for  the  purpose  of 
determining  the  friction  forces)  it  can  be  treated  as  a  cylindrical  pipe 
section.  The  segmental  treatment  of  the  expansion  process  also  accommodates 
the  calculation  of  the  particle  acceleration.  Another  advantage  of  this 
treatment  is  that  it  allows  a  ready  account  of  the  heat  losses  from  the 
wall.  Because  of  the  extremely  high  expansion  temperatures  produced  by 
the  electric  arc  (3200°K),  the  expansion  nozzles  must  be  water-cooled 
with  a  considerable  loss  of  heat  from  the  nozzles. 
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In  principle  the  description  of  a  flow  process  in  terms  of  the 
occurring  fluid  forces  allows  a  completely  consistent  presentation  of  the 
expansion  process.  The  resulting  expansion  Mach  number  is,  together  with 
the  expansion  pressure,  representative  for  the  exit  flow  momentum  and 
constitutes  an  average  flow  Mach  number.  To  obtain  peak  Mach  numbers  an 
assumption  about  the  flow  profile  shape  is  required.  The  present  particle 
acceleration  calculations  use  the  average  Mach  number  given  by  the  flow 
momentum.  Therefore,  the  resulting  particle  velocities  are  also  average 
values . 

The  practical  difficulty  with  applying  the  present  calculation  method 
lies  in  the  choice  of  the  wall  friction  coefficient.  Wall  friction  is  a 
complex  function  of  flow  Reynolds  number,  flow  Mach  number,  and  flow 
history.  The  primary  purpose  of  the  present  calculations  was  to  compare 
particle  velocities  obtained  with  different  nozzle  geometries.  Since 
deviations  from  accurate  friction  values  enter  the  comparison  only  as 
secondary  effects,  constant  wall  friction  coefficients  which  distinguish 
only  between  subsonic  and  supersonic  flow  conditions  were  chosen.  On  the 
same  grounds  the  particle  drag  coefficient  was  assumed  constant,  in  this 
case  for  the  entire  range  of  the  expansion. 

Here  the  wall  friction  is  expressed  in  terms  of  the  pipe  friction 
coefficient.  Since  this  coefficient  is  commonly  referred  to  average 
flow  conditions  derived  from  flow  continuity,  the  Mach  number  used  here 
satisfies  continuity  and  is  particularly  suitable  for  determining  wall 
friction  by  means  of  the  pipe  friction  coefficient. 

For  the  purpose  of  comparing  the  effect  of  nozzle  geometries  on  the 
obtainable  particle  velocities  it  was  not  considered  necessary  to  include 
mutual  interference  between  particles  and  the  influence  of  the  particles 
on  the  gas  flow. 
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SECTION  II 

NOZZLE  FLOW  CALCULATIONS 

1.  FLOW  EXPANSION  IN  TERMS  OF  FLUID  FORCES 

The  balance  between  flow  momentum  and  fluid  pressure  forces  in  a 
duct  with  varying  cross-sectional  area  forms  the  basis  for  the  following 
derivations.  For  a  duct  with  inlet  cross-sectional  area  A-j  and  exit 
cross-sectional  area  A2  as  pictured  in  Figure  1  for  supersonic  expansion 
in  a  parallel  sidewall  nozzle,  the  relation  between  inlet  momentum  I-j 
and  exit  momentum  I2  can  be  written  with  p-j  and  p2  denoting  inlet  and 
exit  pressure 

Pi  ~  Pp 

^  =  l2  -  A2  (P]  -  p2)  -  (A]  -  Az)  -  2  1  (1) 

This  equation  is  equally  applicable  to  a  subsonic  expansion  nozzle  by  a 
consis;ent  use  of  the  signs.  Actually,  the  equation  is  valid  for  any 
type  of  adiabatic  duct  flow  such  as  a  sudden  duct  enlargement  or  a 
supersonic  shock  diffuser.  The  fluid  forces  acting  in  the  axial 
direction  as  a  result  of  the  pressure  difference  in  the  flow  between 
the  duct  walls  and  the  exit  plane  are  expressed  here  with  the  help  of 
a  kind  of  shape  factor  designated  in  Equation  1  with  i.  It  is  introduced 
into  the  calculations  by  the  substitution 


i.  e.,  the  pressure  force  is  calculated  in  a  linearized  way  by  multi¬ 
plying  the  wall  area  projected  in  the  axial  direction  with  the  median 
value  of  the  pressure  difference  existing  between  the  inlet  and  exit 
of  the  duct  element.  The  factor  i  serves  as  a  means  to  account  for  the 
difference  between  the  true  integral  and  the  linearized  calculation. 
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For  a  straight  line  pressure  distribution  in  a  duct  configuration  as  shown 
in  Figure  1,  the  factor  i  is  one.  In  general,  pressure  distributions  in 
expansion  nozzles  are  different  from  a  straight  line  and  vary  in  shape 
with  expansion  conditions  and  also  with  the  shape  of  the  duct  cross  section. 
The  further  development  of  Equation  2  leads  to  an  important  involvement 
of  the  factor  i  in  characterizing  flows  in  ducts  for  the  present  calcula¬ 
tion  method.  Setting  Ag/A-j  =  t.  Equation  1  can  be  written 

1 1  =  1 2  +  (p2  -  Pi  )*2(k + 1  ‘  k  {3) 


The  expression  in  parentheses,  abbreviated  in  the  following  by  writing 


(4) 


becomes  the  sole  criterion  for  characterizing  the  type  of  adiabatic  duct 
flow  to  which  Equation  3  is  applied.  For  instance,  for  a  constant  area 
supersonic  diffuser  it  is  x  5  1  since  t  =  1.  For  a  sudden  enlargment  of 
a  duct  from  A^  to  A?,  it  is  t  =  1/t  as  can  be  derived  from  Equation  3. 
Equation  3  can  also  be  extended  to  account  for  a  second  fluid  mixing 
with  the  first  one  as  it  occurs  in  an  ejector  pump.  For  ejector  mixing 
t  can  again  be  treated  as  a  constant  known  from  experiments.  Supersonic 
diffusers  and  ejectors  have  been  treated  in  connection  with  other 
investigations  (Reference  1  and  2).  For  an  adiabatic  expansion  process 
(which  is  of  interest  here)  r  must,  in  general,  be  determined  in  each 
case  separately  from  the  expansion  conditions  (Reference  2).  However, 
as  pointed  out  in  the  Introduction,  for  the  present  case  the  calculations 
are  applied  in  a  stepwise  fashion  to  small  segments  of  the  nozzle  and  x 
can  also  be  assigned  a  unique  value.  If  the  considered  nozzle  segment  is 
short  enough,  in  particular  if  the  pressure  drop  across  the  segment 
is  small,  the  pressure  distribution  over  the  length  of  the  segment 
approaches  a  straight  line,  i.  e.,  i  approaches  unity  with  the  influence 
of  the  nozzle  cross-section  shape  in  i  also  disappearing.  Then 
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At  the  end  of  this  Section  it  will  be  shown  that  this  simplified  approach 
leads  to  very  accurate  results  for  the  overall  expansion  in  a  nozzle  if 
applied  successively  to  small  segments  of  the  expansion  nozzle.  Equation 
3  will  be  further  developed  with  t  considered  as  a  known  magnitude 
according  to  Equation  5  by  the  geometry  of  the  contemplated  segment 
of  the  nozzle. 

2.  INCLUSION  OF  THE  WALL  FRICTION  FORCES 

Considering  the  generally  known  relation 

I  =  A  p  y  M2  (6) 

with  y  as  the  ratio  of  specific  heats  and  M  as  the  Mach  number  of  the  gas 
flow  Equation  3  becomes  with  the  inclusion  of  the  wall  friction  forces 

Ff 

A1  P'1  yMi  =  A2‘P2'Y‘M2  +  T  A2  (p2  "  pl^  +  Ff 


The  wall  friction  forces  can  be  expressed  in  terms  of  the  common  pipe 
friction  coefficient  in  the  following  form 


or  introducing  the  flow  Mach  number 

cf  L  2 

F^  -  2  o  p2  T  ^2  A2  (9) 

For  the  pipe  as  a  constant  geometry  duct  the  inlet  is  the  obvious 
choice  as  the  reference  location  for  the  friction  coefficient  since  the 
flow  conditions  are  normally  given  for  this  location  in  case  of  a  normal 
pipe.  For  a  duct  with  changing  cross  section  as  treated  here,  the  choice 
of  reference  cross  section  depends  on  the  significance  of  the  flow  con¬ 
ditions  for  the  overall  wall  friction  in  the  chosen  duct  segment.  The 
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present  calculation  method  has  been  originally  derived  for  treating 
ejector  systems.  For  the  ejector,  the  mixed  flow,  i.e.,  the  exit  flow  is 
more  influential  for  the  overall  wall  friction  than  the  inlet  flow  where 
only  the  slow  moving  induced  flow  is  in  touch  with  the  wall.  The  exit 
of  the  ejector  mixing  section  is  therefore  the  logical  choice  for  the 
friction  coefficient  reference  location  (Reference  1).  For  an  expansion 
nozzle  the  choice  becomes  more  arbitrary,  particularly,  if  the  calcula¬ 
tions  are  applied  to  very  short  segments  of  the  nozzle  as  proposed  here. 
To  allow  a  consistent  use  of  the  originally  derived  equations  the  exit 
of  the  segment  was  chosen  as  the  reference  location. 


Equation  7  can  now  be  transformed  into 


y-M 


M2,CfL 
Y  ^2  2D 


+  1  )  +  T 


=  B 


(10) 


This  equation  contains  two  unknowns:  the  exit  Mach  number  and  the 
pressure  ratio  across  the  expanding  flow.  Conservation  of  mass  and 
energy  as  expressed  by  the  following  relations 


(11) 


p  = 


P_ 

RT 


(12) 


(13) 


allows  one  to  write  an  additional  condition  for  determining  the  two  unknowns. 
Equations  11  through  13  yield  the  following  relation 


h2(i  + 


Mi  ^ 
t  P2 


(14) 
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Solving  Equation  10  for  /p^  and  introducing  it  into  Equation  14  leads 


to  the  basic  equation  of  our  flow  problem 


.1  ?  V2 

( 1  +  Mg) 

o  c  fL 

Y  M2(g^-  +  1)  +  x 


! 

! 


(is)  : 

i 


This  equation  can  be  solved  for  the  exit  Mach  number 


M 


2 

2 


1 


.f2,2£  +  l.n*ot] 

Cay  -  Y  +  1 


(16) 


using  the  abbreviations 


V.  ~l_ 

c  =  1  +  2D~ 

(17) 

a  =  2  yE2C 

(18) 

The  magnitude  E  in  Equation  18  stands  for  the  expression  on  the  right 
side  of  Equation  15.  The  expansion  pressure  ratio  with  Equation  10 
becomes 


P2 


Pi 


+ 


c  + 


T 


T 


(19) 


Equations  16  through  19  are  applicable  to  any  adiabatic  duct  flow 
including  flow  mixing,  in  which  case  only  the  right  side  of  Equation  15 
assumes  a  different  form.  The  positive  root  in  Equation  16  yields  the 
subsonic  and  if  applicable,  the  negative  root  the  supersonic  solution 
of  the  flow  problem.  The  pipe  friction  coefficient  in  Equation  17  is  a 
commonly  known  function  of  the  flow  Reynolds  number  and  the  wall  surface 
roughness.  As  indicated  before  x  is  found  from  Equation  5. 
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Equation  16  can  be  applied  to  any  two  cross  sections  of  a  complete 
supersonic  expansion  nozzle  consisting  of  a  subsonic  and  a  supersonic 
portion.  The  determination  of  t  requires  a  more  complex  procedure  if  the 
restriction  to  small  segments,  i.  e.,  to  small  flow  changes  is  not 
maintained.  The  necessary  relations,  which  are  of  no  interest  here, 
have  been  derived  previously  in  connection  with  other  investigations 
(Reference  2). 

In  Figure  2,  Equations  16  and  19  are  evaluated  for  zero  wall  friction 
losses  using  Equation  5  for  determining  t.  Slightly  different  curves 
result  depending  on  the  choice  of  the  segment  length  aL.  For  comparison, 
results  obtained  with  exact  t  values  are  also  entered.  As  expected  on 
the  basis  of  the  definition  of  t,  the  exact  results  are  independent  of 
the  segment  length.  They  are  also  in  exact  agreement  with  the  values 
obtained  from  the  common  isentropic  expansion  relations.  From  Figure  2 
one  sees  that  for  practical  purposes  the  segmental  calculation  method 
using  a  simplified  -r-value  yields  very  satisfactory  results  for  a  proper 
choice  of  the  segment  lengths. 
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SECTION  III 

THE  PARTICLE  ACCELERATION  PROCESS 

1.  SEGMENTAL  CALCULATION  PROCESS 

The  acceleration  of  a  particle  in  a  gas  flow  can  be  based  on  the 
consideration  that  the  work  done  on  the  particle  by  the  drag  forces 
over  a  certain  path  length  must  appear  again  as  kinetic  energy  addition 
to  the  particle.  Using  the  notations  given  in  Figure  3a  this  condition 
can  be  expressed  by  the  following  relation 

(v-j  -  if  Cq  (^-)^  ir  aL  =  Pp  j  Tt  (^-)  [  (vp  -  (vp  ]  (20) 

or  with  some  rearrangement 

4L  ('i  -  ,2>2  ■  [  (»2>2  -  <'2>2  ]  <21> 

In  this  approach  any  change  of  the  gas  velocity  within  the  particle  path 
AL  is  neglected.  It  will  be  shown  later  that  this  simplified  approach 
is  still  admissible  for  the  case  of  the  expanding  flow,  i.e.,  for  an 
increasing  gas  velocity,  if  AL  is  chosen  sufficiently  small. 

Equation  21  can  be  developed  further  by  solving  it  for  v£  and 
dividing  both  sides  with  v!j 


The  velocity  ratio  obtained  with  this  relation  is  the  exit  velocity  ratio 
of  the  contemplated  segment  and  thus  the  inlet  velocity  ratio  for  the 
following  segment. 
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To  check  the  accuracy  of  Equation  22,  a  more  elaborate  relation 
which  takes  the  change  of  the  flow  velocity  in  the  segment  into  account, 
was  also  derived.  The  pertinent  acceleration  model  is  shown  in  Figure  3b. 
The  conditions  for  this  model  can  be  written  in  the  following  differential 
form 

(v-j  -  v2)2  c p  2^  (f)2  *  ds  =  \  Pp  J  TT  (|)3  [(v2  +  dv2)2  -  v2j  (23) 
or  solved  for  ds 
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it  can  be  neglected  and  the  equation  becomes  with  dividing  numerator  and 
2 

denominator  by  (vp 


ds 


v?  v? 

2  k  4  d(4) 

v-j  vp 


'  ■  4 +  + 2<^>s  * 


(29) 


setting  s  =  x 

v2/vl  =  * 


the  following  final  differential  equation  can  be  formed 


= 

dx 


1  [4- 

k  L2y 


-  1  + 


1  +  *  (L.)  +  x 

2  y  VALJ 


_  (^-)2  .  x(— )] 

2y  1aL;  MaLJJ 


(30) 


This  equation  can  be  solved  by  readily  available  computer  methods  (For  the 
present  evaluations  the  Runge-Kutta  method  on  the  Hewlett-Packard  Calcula¬ 
tor  9820a  was  used).  Figure  4  compares  the  simplified  method  given  by 
Equation  22  with  the  more  accurate  one  expressed  in  Equation  30  for 
the  same  segment  length  in  each  case.  One  sees  that  the  simplified  method 
gives  practically  the  same  results  as  the  more  refined  calculations.  The 
agreement  begins  to  suffer  with  the  particle  speed  approaching  the  gas 
speed  as  can  be  seen  from  the  upper  curve  with  the  smaller  particle  to  gas 
density  ratio,  where  the  particle  speed  reaches  about  75%  of  the  gas  speed. 
Since  for  the  present  calculations  (Section  V)  the  velocity  ratios  stay, 
in  general,  below  this  value,  the  simplified  method  was  exclusively  used 
in  these  calculations. 

2.  SEGMENT  SELECTION 

A  simple  relation  is  derived  which  allows  a  convenient  selection 
of  the  segment  length  for  the  stepwise  calculation  of  the  expansion 
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process  in  a  nozzle  with  a  circular  cross  section.  From  the  schematic 
drawing  of  a  portion  of  an  axisymmetric  supersonic  expansion  nozzle 
shown  in  Figure  5  one  can  derive  the  following  relations 

D,  2  A, 

(—1  =  —  I 

u2  h2 


0^  =  -  2  L  tan  B 


These  two  relations  can  be  combined  to  yield 


2 

(1  -  2  t  tan  6) 
u2 


In  this  way  the  segment  cross  section  ratio,  as  it  enters  the  nozzle  flow 
equations,  is  related  to  readily  given  nozzle  geometry.  For  the  subsonic 
portion  of  an  expansion  nozzle  the  corresponding  equation  is 


(1  +  2  tan  B) 
u2 


12 


AFWAL-TR-81  -31 05 


SECTION  IV 
NOZZLE  HEAT  LOSSES 


The  stepwise  calculation  of  the  gas  expansion  process  provides 
ready  conditions  for  the  inclusion  of  heat  losses  through  the  nozzle 
wall.  For  the  heat  loss  from  a  small  nozzle  segment  one  can  write 


AQLoss  =  h-D  '*  AL-At-AT 


(35) 


This  equation  will  be  treated  in  more  details.  The  heat  transfer 
coefficient  h  in  this  relation  is  conveniently  obtained  with  the  help 
of  the  Reynolds  analogy,  which  states  (Reference  3) 


Nu 


Re 


^f^wall 

2 


(36) 


Then  the  heat  t  ansfer  coefficient  follows  using  the  well  known 
relation 

(c  )  - 

cf  wall  '  4  (37) 


h 


k  v 


llflpi  pe 
8 


(38) 


Thus,  we  can  introduce  our  previously  used  pipe  friction  coefficient. 

The  time  element  At  in  Equation  35  stands  for  the  residence  time  of  the 
gas  in  the  segment  and  is  given  by 

At  =  (39) 


The  temperature  difference  aT  in  Equation  35  is  given  by  the  relation 

ST  ■  'VoW-  <4t» 

where  «rec  is  the  wall  temperature  recovery  factor.  The  remaining  factors 
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in  Equation  35  are  the  segment  Diameter  D  and  the  segment  length  aL. 
Therefore,  Equation  35  using  Equations  38  to  40  becomes 


AQ 


loss 


Du  ( aL) 


2 


(V 


rec 


-  T 


wal  1 


(41) 


Due  to  the  heat  loss  given  by  this  equation  the  total  gas  temperature 
drops  within  the  nozzle  length  aL 


<T2>0 


=  (T.) 

1  0 


AO, 

1  OSS 

c  •  m 
P 


(42) 


In  this  equation  m  is  the  gas  mass  passing  through  the  segment  within  the 
time  element  At.  For  the  segment  exit  temperature  one  can  write  more 
expl icitly 


<T2>  o 


<T1>  o 


k  c , 


"V 


r  „  W  -  TwaiP 


(43) 


The  nozzle  wall  temperature  in  this  equation  must  still  be  determined. 
Assuming  water-cooled  walls,  i.  e.,  a  constant  temperature  on  the  outer 
nozzle  wall,  one  can  write  based  on  the  common  heat  conduction  relation 


wall 


Tlux 


Xu 


water 


(44) 


The  heat  flux  occurring  in  this  equation  is  readily  obtained  from  the 
previously  derived  relation  for  the  heat  loss  from  the  segment  by 
writing 


AQ 


Loss 


'flux 


S-At 


(45) 


Since  AQ^qss  is  a  function  of  the  wall  temperature,  Q  is  also  a 

function  of  the  wall  temperature.  To  overcome  this  functional  difficulty 
the  nozzle  calculations  are  started  with  an  estimated  wall  temperature 
for  the  first  segment  of  the  nozzle. 
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Applying  Equation  43  to  consecutive  segments  of  the  nozzle  the 
drop  in  total  temperature  over  the  whole  nozzle  length  can  be  determined. 
The  correction  of  the  gas  density  as  it  is  needed  for  calculating  the 
particle  acceleration  can  now  be  found  from  the  total  temperature 
witi  the  pressure  in  the  flow  and  the  flow  Mach  number  known  from  the 
previously  derived  flow  calculations  (Equations  16  and  19). 
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SECTION  V 

CALCULATION  RESULTS 

Figures  6  through  8  show  the  flow  conditions,  the  particle  velocities, 
and  the  heat  losses  calculated  for  a  number  of  axisymmetric  nozzle  con¬ 
figurations  for  one  set  of  operating  conditions,  which  is  only  varied 
for  Figure  8  where  the  particle  injection  point  is  moved  downstream 
of  the  nozzle  throat.  The  geometry  given  by  Figure  6  for  the  subsonic 
and  by  Figure  7a  for  the  supersonic  nozzle  portion  is  that  of  an  exper¬ 
imental  acceleration  nozzle.  All  configurations  have  the  same  subsonic 
nozzle  geometry  as  given  in  Figure  6,  which  provides  the  calculation 
results  for  the  subsonic  conditions.  Figures  7a  to  c  give  the  results 
for  the  supersonic  nozzle  portion  for  various  nozzle  geometries.  The 
geometries  are  identified  in  the  figures  by  a  plot  of  the  nozzle  diameter 
over  the  nozzle  length,  which  serves  as  abscissa  in  all  presentations 
of  the  calculation  results.  The  nozzle  wall  angle  is  also  listed  in 
the  figures. 

Nozzle  wall  friction  and  particle  drag  coefficients  used  for  the 
present  calculations  are  estimates  based  on  generally  available  informa¬ 
tion  on  the  subject.  Typical  Reynolds  numbers  for  the  present  type  of 
flow  conditions  are  listed  in  Table  1. 


TABLE  1 

REYNOLDS  NUMBERS  FOR  FLOW  CONDITIONS 


Nozzl e  Pos  .  L 

Re 

P 

Ref 

14  mm 

940 

104000 

20 

1300 

117000 

28 

870 

77000 

41 

700 

68000 

65 

530 

59000 

116 

350 

49000 
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TABLE  1  (Continued) 

REYNOLDS  NUMBERS  FOR  FLOW  CONDITIONS 


Nozzle  Pos.  L 

% 

Ref 

241 

180 

37000 

399 

100 

29000 

602 

60 

22000 

806 

27 

14000 

For  the  subsonic  portion  of  the  expansion  nozzles  a  friction 
coefficient  of  0.02  (pipe  friction  coefficient)  was  chosen.  For  the 
supersonic  portion  the  corresponding  value  is  0.015.  Nozzle  flow 
calculations  presented  in  Reference  4  provided  some  special  guidelines 
for  these  choices.  An  analysis  on  wall  friction  losses  given  in 
Reference  2  also  influenced  the  choices.  Though  experiments  with 
expansion  nozzles  having  pertinent  geometries  were  carried  out,  no 
direct  comparison  of  actual  gas  or  particle  speeds  with  calculated 
values  were  possible  due  to  the  difficulties  in  reliably  measuring 
velocities  under  the  applicable  extreme  test  conditions.  However,  in 
the  present  analysis  the  Reynolds  analogy  between  friction  and  heat 
losses  was  used  to  determine  nozzle  heat  losses  and  subsequently  gas  exit 
temperatures.  Thus,  a  comparison  of  calculated  and  experimental  gas 
temperatures  provides  a  clue  for  the  validity  of  the  calculations. 
Circumstances  did  not  allow  exact  gas  temperature  measurements  in  the 
present  tests.  However,  from  the  glow  color  of  an  uncooled,  190  mm  long, 
stainless  steel,  guide  tube  attached  to  the  nozzle  exit  an  estimate  of 
the  gas  temperature  was  possible.  The  yellow-white  glow  color  of  the 
exit  portion  of  the  guide  tube  suggested  a  tube  color  temperature  of 
1400°  to  1500°K.  The  calculated  gas  temperature  at  the  nozzle  exit 
(L  =  660  mm  in  Figure  7a)  is  1600°K.  Since  the  heat  losses  from  the 
tube  (:  0.06  KW/cm2 )  were  small  compared  to  the  heat  flux  from  the 
gas  (:  0.7  KW/cm2),  the  tube  end  which  is  fairly  far  removed  from  the 
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tube  attachment  point,  can  be  assumed  to  have  nearly  reached  gas  recovery 
temperature.  Therefore,  calculation  and  experiment  agree  within  the  order 
of  I00°K,  a  value  which  allows  one  to  conclude  that  the  nozzle  wall 
friction  assumptions  are  reasonable. 

For  the  part  of  the  expansion  nozzle  where  most  of  the  particle 
acceleration  occurs,  the  particle  Reynolds  number  ranges  from  about  1300 
near  the  throat  down  to  about  100  at  L  =  400  mm.  The  drag  coefficient  for 
a  sphere  in  this  regime  ranges  from  about  0.4  to  1.0.  Since  the  actually 
used  particles  {graphite  pellets)  deviated  from  a  spherical  shape  and  had 
no  smooth  surface,  a  drag  coefficient  near  the  upper  limit  was  chosen 
with  a  value  of  0.8.  For  simplicity  this  value  was  also  applied  to  the 
nozzle  exit  region  where  the  particle  approaches  the  Stokes  flow  regime. 

Figures  7a  to  c  show  that  the  obtainable  particle  velocity  is  not 
greatly  affected  by  the  nozzle  geometries  considered  in  these  figures. 

The  particle  to  gas  velocity  ratio  never  exceeds  a  value  of  around  0.5 
except  in  the  regions  where  the  gas  velocity  drops  from  its  peak  value. 

In  all  cases  the  particle  velocity  increases  very  slowly  in  most  of  the 
downstream  portion  of  the  supersonic  expansion  region. 

For  selecting  the  best  nozzle  geometry  the  obtainable  particle  veloci¬ 
ties  must  be  compared  for  the  expansion  pressure  ratio  given  by  the 
experimental  conditions.  In  the  present  case  this  ratio  is  determined 
by  the  permissible  arc  heater  operating  pressure  (100  kg/cm^)  and  ambient 
conditions  (1  kg/cm  )  preferred  for  carrying  out  the  erosion  experiments. 
Resulting  particle  velocities  are  listed  in  Table  2.  This  table  gives 
an  overview  for  the  general  possibilities  of  accelerating  80-micron 
particles  in  an  expanding  gas.  Cases  No.  1  to  3  pertain  to  the  conditions 
given  in  Figures  7 a  to  c .  For  the  largest  pressure  ratio  possible  with 
the  present  experiments(p/p0  =  0.01),  the  maximum  particle  velocity  is 
obtained  in  Case  No.  2  with  a  0.5°  wall  angle  for  the  supersonic  nozzle 
portion.  This  velocity  of  2002  m/sec  is  only  slightly  better  than  the 
ones  obtained  in  Cases  No.  1  and  3.  The  particle  velocities  given  for 
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pressure  ratio  0.02  are  indicative  for  the  conditions  where  the  nozzles 
are  cut  off  before  ambient  pressure  is  reached.  The  particle  velocities 
are  only  slightly  reduced  in  these  cases. 

The  remaining  cases  in  Table  2  serve  to  demonstrate  the  practicality 
of  choosing  the  small  size  nozzles  given  by  Cases  No.  1  to  3  and  dealt 
with  in  details  in  Figures  7a  to  c.  Cases  No.  4  to  7  show  that  higher 
particle  velocities  can  be  obtained  by  increasing  the  nozzle  size  and 
consequently  the  nozzle  mass  flow.  For  the  same  particle  residence  time 
in  the  gas  the  relative  wall  friction  influence  is  reduced.  By  doubling 
the  nozzle  cross-sectional  areas  (Cases  No.  4  and  5)  some  gain  in  particle 
velocity  is  obtained.  In  view  of  the  associated  increase  in  cost  of 
already  expensive  equipment  and  operating  medium  the  worthwhileness  of 
the  improvement  is  questionable.  A  further  doubling  of  the  nozzle  cross- 
sectional  areas  (Cases  No.  6  and  7)  produces  a  further  increase  in 
particle  velocity.  However,  cost  may  become  prohibitive  in  these  cases. 

Cases  No.  8  and  9  show  examples  with  air  as  the  operating  medium. 

Due  to  its  higher  molecular  weight  air  cannot  be  expanded  to  as  high  a 
velocity  as  helium,  however,  its  higher  density  provides  a  better  drag 
coupling  between  particle  and  gas  (Equation  22).  Thus,  air  is  of  some 
interest  for  accelerating  80-micron  particles.  Case  No.  8  is  the  same 
as  Case  No.  3  with  air  replacing  helium  as  operating  medium.  Only  the 
air  velocities  are  listed  for  air.  Survey  calculations  indicated  that 
due  to  the  superior  coupling  in  this  case  the  particle  velocities 
approach  closely  the  air  velocities.  This  closeness  makes  the  present 
method  of  calculating  particle  velocities  invalid.  The  air  velocities 
listed  in  Table  2  serve  as  upper  limit  for  the  particle  velocities. 

Even  the  air  velocities  stay  substantially  below  the  particle  velocities 
obtained  with  helium. 

Case  No.  9  deals  with  the  practical  possibility  of  using  a  large  scale 

p 

air  expansion  nozzle  with  increased  plenum  pressure  (200  kg/cm  ).  To  stay 
near  practical  limits  the  plenum  temperature  is  reduced  to  2000°K.  Com¬ 
pared  to  Case  No.  8  the  air  velocities  are  substantially  increased  but 
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still  appreciably  below  the  particle  velocities  obtained  with  the  much 
smaller  helium  operated  nozzles. 

The  drop  in  gas  velocity  in  the  downstream  portion  of  the  supersonic 
nozzle  is  caused  by  the  large  drop  in  the  gas  stagnation  temperature  due  to 
the  heat  loss  from  the  nozzle.  The  figures  show  that  the  expansion  Mach 
number  still  increases  throughout  the  nozzle  though  mostly  at  a  very  small 
rate  caused  by  the  accumulative  effect  of  wall  friction.  An  increasing 
Mach  number  is  a  requirement  for  maintaining  supersonic  flow  stability 
in  the  nozzle. 


In  the  throat  region  the  heat  flux  to  the  nozzle  wall  reaches  a 
2 

maximum  of  about  10  KW/cm  .  This  enormous  heat  flux  calls  for  particular 
attention  in  the  design  of  the  nozzle  cooling  system.  The  step  down  in 
the  heat  flux  seen  in  the  curves  at  the  nozzle  throat  is  simply  caused  by 
the  change  of  the  wall  friction  coefficient  from  0.02  for  the  subsonic 
nozzle  portion  to  0.015  for  the  supersonic  portion. 

Figure  8  shows  the  case  given  in  Figure  7a  with  the  point  of  particle 
injection  moved  downstream  of  the  throat.  In  spite  of  the  delayed  injection 
the  final  particle  speed  is  reduced  only  slightly  from  the  case  shown  in  Fiqure 
7a.  As  a  general  observation  one  can  say  that  in  all  cases  treated  here, 
the  particles  are  very  rapidly  accelerated  to  about  half  the  gas  velocity 
but  further  acceleration  then  becomes  increasingly  sluggish. 
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SECTION  VI 

FLOW  PROFILE  CONSIDERATIONS 


1.  GENERAL  DISCUSSION 

The  present  flow  calculation  method  is  set  up  to  find  for  a  given 
duct  the  flow  momentum  at  the  exit  of  the  duct  from  the  flow  conditions 
at  the  inlet.  The  result  is  expressed  in  the  form  of  an  effective  Mach 
number.  This  Mach  number  represents  average  flow  conditions  which  in 
accordance  with  its  derivation  fulfills  the  condition  of  continuity  in 
the  flow  through  the  duct.  Since  the  definition  of  the  pipe  friction 
coefficient  is  based  on  the  use  of  such  average  flow  conditions  as  a 
reference,  the  effective  Mach  number  is  the  right  choice  for  determining 
the  friction  forces  (Reference  3). 

For  the  particle  acceleration  it  is  of  interest  to  obtain  infor¬ 
mation  about  the  particle  velocity  distribution  over  the  flow  cross 
section,  in  particular,  about  the  peak  velocities  of  the  particles. 

In  a  first  approximation  the  particle  velocity  distribution  at  the  nozzle 
exit  can  be  considered  as  directly  proportional  to  flow  velocity  dis¬ 
tribution,  i.e.,  it  is  sufficient  to  determine  the  particle  velocity 
distribution  from  the  flow  velocity  distribution  by  simply  using  the 
particle  to  gas  velocity  ratio  found  in  the  calculations  for  the  nozzle 
exit  as  a  proportionality  factor  for  the  two  distributions.  For  finding 
the  peak  particle  velocity  by  this  process  a  method  is  presented  in  the 
following  which  allows  one  to  derive  the  peak  Mach  number  from  the 
effective  Mach  number  for  a  given  flow  profile  shape. 

2.  PEAK  MACH  NUMBER 

If  the  flow  profile  shape  is  given,  the  peak  Mach  number  can  be 
uniquely  related  to  the  effective  Mach  number  by  the  requirement  that 
for  each  Mach  number  distribution  the  flow  momentum  must  be  the  same 
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where  y  and  p  are  considered  constants.  The  flow  profile  shape  is  con¬ 
sidered  to  be  given  in  the  form  of  a  Mach  number  distribution  over  the  radius 
of  a  circular  cross  section.  Determining  the  peak  Mach  number  is  a  simple 
geometric  problem.  For  the  present  purpose  the  shape  of  the  profile  is 
presented  in  terms  of  the  square  of  the  Mach  number  and  assumed  to  be 
the  frustrum  of  a  cone.  Then  the  specific  problem  is,  as  illustrated 
in  Figure  9,  to  find  the  height  of  the  frustrum  which  has  the  same 
volume  as  a  cylinder  the  height  ( hs+h^ )  of  which  is  given.  Practical 
relations  are  derived  in  the  following. 

2 

The  M  -  profile  is  given  by  the  two  ratios 


For  establishing  a  relation  between  h$  and  h 
is  selected 


an  auxiliary  parameter 

(48) 


The  following  geometric  conditions  must  be  fulfilled 

1.  from  Equation  46  Volume  V-|  =  Volume 

2.  from  similarity  ^o  ^o  +  ^1  ^o  +  ^1  +  ^2  (49) 

do  '  dl  d2 
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From  these  conditions  one  can  derive 


£>3 

al 


•  3<df> 


o  *  >  ) 


The  solution  to  this  relation  is 


d,  i  ,3 

jj-  *  2  cos [■j  arccos( - p5-)] 


With  Equation  51  the  auxiliary  parameter  a  is  now  directly  connected 
to  the  above  given  M  -  profile  description  considering  that 


Simple  geometric  similarity  now  allows  all  magnitudes  of  the  profile 
scheme  in  Figure  9  to  relate  to  each  other.  In  particular,  Equation  47 
can  now  be  conveniently  written 


Meff  2 
max 


+  1 


h.  h 

_L  +  i  +  _L 

h  h2 


where  the  occurring  ratio  values  on  the  right  side  are  related  to  the  given 
profile  conditions  by 


1  -  a 

i 

dl 


hs  hl 
ht  h2 


+  1 


1  - 
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Equations  51  through  55  allow  one  to  prepare  a  plot  as  shown  in 
Figure  10.  The  flow  profiles  shown  in  this  figure  give  the  direct  Mach 
number  profiles  over  the  duct  diameter  instead  of  the  square  values  as 
originally  expressed  in  Figure  9.  Figure  10  covers  all  geometric  condi 
tions  with  which  flow  profiles  can  be  expressed.  Actually  occurring 
flow  profiles  can  be  identified  with  those  occurring  along  very  low 
values  of  M,  .,,/M  except  for  high  d./d,  values,  where  large  Mach 

number  ratios  are  possible  characterizing  the  flow  profile  start-up  at 
the  inlet  to  a  duct.  Note:  in  Figure  10  notation  dQ  is  changed  to  d-j . 
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SECTION  VII 
CONCLUSIONS 

The  acceleration  of  80p  diameter  particles  in  expanding  helium  to 
near  60%  gas  velocity  for  the  purpose  of  re-entry  erosion  testing  requires 
the  use  of  expansion  nozzles  which  have  a  length  in  the  order  of  100 
nozzle  throat  diameters.  Wall  friction  incurred  in  such  nozzles  becomes 
a  limiting  factor  for  the  effectiveness  of  this  acceleration  process. 

To  account  for  the  friction  influence  in  the  design  of  these  nozzles  the 
expansion  process  is  conveniently  expressed  in  terms  of  the  occurring 
fluid  forces.  This  procedure  allows  a  ready  inclusion  of  the  friction 
forces  which  by  themselves  can  be  determined  with  the  help  of  known  wall 
friction  coefficients.  This  calculation  method  yields  a  particularly 
accurate  account  cf  the  friction  effects  if  applied  consecutively  along 
the  nozzle  to  small  nozzle  segments.  This  stepwise  calculation  procedure 
also  accommodates  the  determination  of  the  heat  losses  from  the  nozzle 
and  provides  a  basis  for  calculating  the  particle  acceleration  process. 

The  segment  length  can  be  adjusted  to  the  rapidity  of  the  flow  changes 
as  it  is  particularly  necessary  in  the  nozzle  throat  region  due  to  the 
extreme  sensitivity  of  the  flow  to  cross  section  changes  in  this  region. 
Very  small  segments  must  be  chosen  to  stay  within  the  validity  of  the 
assumptions  made  for  the  calculations. 

The  present  flow  calculation  process  does  not  include  detail  boundary 
layer  considerations.  Instead  it  accounts  for  viscous  effects  in  terms 
of  a  wall  friction  coefficient.  In  spite  of  this  simplification  a 
consistent  description  of  the  flow  process  is  possible  as  long  as  the 
chosen  friction  coefficient  is  sufficiently  representative  for  the 
viscous  conditions  in  the  flow  or  if  the  coefficient  is  only  used  as  a 
parameter  to  demonstrate  its  influence  on  the  flow.  For  the  present 
purpose  of  comparing  particle  velocities  obtained  with  only  slightly 
changed  nozzle  geometries  the  exact  representation  of  the  viscous 
conditions  was  not  considered  essential  and  constant  friction  coefficient 
values  only  different  for  subsonic  and  supersonic  flow  were  chosen. 
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The  end  results  of  the  present  flow  calculations  give  the  flow  exit 
momentum  expressed  in  terms  of  an  average  flow  Mach  number  which  may  also 
be  called  effective  Mach  number.  Assumptions  about  the  flow  profile  are 
necessary  to  determine  the  peak  flow  Mach  number.  For  the  present 
determination  of  particle  velocities  this  effective  Mach  number  was  used. 

In  calculating  nozzle  designs  for  accelerating  small  particles 
(0.08  mm)  to  high  velocities  (2000  m/sec)  one  finds  that  it  is  difficult 
to  obtain  particle  velocities  much  larger  than  50%  of  the  occurring 
maximum  gas  velocities.  Due  to  the  wall  friction  influence  in  the  long 
expansion  nozzles  (0.7°  -  1.4°  opening  angle  for  the  supersonic  nozzle 
portion)  less  than  50%  of  the  isentropic  expansion  velocity  can  be 
obtained  for  the  average  accelerating  gas  velocity.  Since  extreme  high 
gas  temperatures  (arc  heated  gas)  must  be  employed  to  obtain  sufficiently 
high  gas  velocities  (at  an  expansion  pressure  ratio  of  100  and  helium 
as  operating  medium)  the  expansion  nozzles  are  exposed  to  extreme  heat 
loads  (lOKW/cm  in  the  throat  region)  requiring  intense  water  cooling 
of  the  nozzles.  Enthalpy  losses  ranging  from  30%  to  50%  of  the  initial 
enthalpy  occur  during  expansion. 
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Figure  2.  Comparison  of  Segmental  Flow  Calculations  with  Isentropic 

Flow  Conditions  for  Zero  Wall  Friction  (Y  =  1.4).  (For 

clarity,  curve  points  for  the  smaller  segment  length  have  been 
largely  omitted  in  the  supersonic  region.  Differentiating  markings 
of  curve  points  have  been  applied  only  in  the  supersonic  region 
and  only  for  the  Mach  number  curve  where  results  deviate  sufficiently) 
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Flow  Conditions,  Heat  Loss,  and  Particle  Acceleration  for  the 
Subsonic  Portion  of  the  Expansion  Nozzles  Considered  in 
Figures  7a  to  c  for  Their  Supersonic  Nozzle  Portions 


Flow  Conditions,  Heat  Loss,  and  Particle  Acceleration  for 
Various  Supersonic  Nozzle  Configurations  (For  Subsonic 
Nozzle  Portion  See  Figure  6) 
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of  Figure  7 a  with  Particles  Injected  Downstream  of  the 
e  Throat. 


Figure  9.  Schematic  Flow  Profile 
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for  Determining  Peak  Mach  Numbers  from  Effective 
Numbers  for  Given  Flow  Profile  Shapes 


